Light emission from a scanning tunneling microscope: Fully retarded calculation 
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The light emission rate from a scanning tunneling microscope (STM) scanning a noble metal 
surface is calculated taking retardation effects into account. As in our previous, non-retarded theory 
[Johansson, Monreal, and Apell, Phys. Rev. B 42, 9210 (1990)], the STM tip is modeled by a sphere, 
and the dielectric properties of tip and sample are described by experimentally measured dielectric 
functions. The calculations are based on exact diffraction theory through the vector equivalent 
of the Kirchoff integral. The present results are qualitatively similar to those of the non-retarded 
calculations. The light emission spectra have pronounced resonance peaks due to the formation of a 
tip-induced plasmon mode localized to the cavity between the tip and the sample. At a quantitative 
level, the effects of retardation are rather small as long as the sample material is Au or Cu, and 
the tip consists of W or Ir. However, for Ag samples, in which the resistive losses are smaller, 
the inclusion of retardation effects in the calculation leads to larger changes: the resonance energy 
decreases by 0.2-0.3 eV, and the resonance broadens. These changes improve the agreement with 
experiment. For a Ag sample and an Ir tip, the quantum efficiency is « 10 -4 emitted photons in the 
visible frequency range per tunneling electron. A study of the energy dissipation into the tip and 
sample shows that in total about 1 % of the electrons undergo inelastic processes while tunneling. 



PACS numbers: 61.16.Gh, 41.20.Bt, 73.20.Mf, 78.20.Bh 



I. INTRODUCTION 

The scanning tunneling microscope (STM) has over 
the last 15 years developed into a standard instrument 
in surface science. A relatively small, but very interesting 
part of this development is concerned with light emission 
from STM's. This phenomenon, in which the tunneling 
electron transfers energy to a photon p has been observed 
on noble^rnetal surfaces by Gimzewski, Berndt,|-apd co- 
workers,B~El and later by several other groups-Lliiil J-"igh± 
emission was also observed from .semiconductors, E5t£l 
from surfaces covered by molecules,Ej and from magnetic 
surfaces.E3 Other experiments have studied closely re- 
lated phenomena. To mention a few examples, the STM 
tip can be used_as a local detector of photoelectrons in 
photoemission,E2 and if the STM is illuminated by laser 
light, the rectification current can be large enough to 
generate STM topographs.!!!] Moreover, the interaction 
between propagating surface plasmons and an STM tip 
has been investigated.tEl 

Light emission from noble-metal surfaces-i-as inspired 
a certain amount of theoretical activityEJcIl The most 
striking experimental results, the high (in this context) 
quantum efficiency of the process (up to 10 -3 photons 
per electron) and the characteristic resonances in the 
light emission spectra, can be understood in the follow- 
ing way:E!Tca When the tip-sample separation is no more 
than 5-10 A, the surface plasmons on the two surfaces 
interact quite strongly and form an interface-plasmon 
mode. This is illustrated schematically in Fig. |l|. The 
resonance in the light emission spectrum occurs at the 
frequency for which half a wavelength of the interface 
mode fits into the "cavity" between tip and sample. The 
charge oscillations associated with the interface plasmon 



have opposite signs on the two electrodes. This has two 
major consequences: (i) The charge on one surface gives 
a field that polarizes the second surface so that the polar- 
izing field at the first surface increases, etc. Thus, given 
the charge configuration of the interface mode, the en- 
hancement of the electromagnetic vacuum fluctuations is 
a quite natural consequence, (ii) Because the charges on 
the tip and the sample are opposite, the interface plas- 
mon is locally charge-neutral. This weakens the restor- 
ing forces compared with the case of a surface plasmon 
at an isolated interface, and the resonance peak is red- 
shifted relative to the surface plasma frequency. From 
this argument follows also a geometry-dependence of the 
resonance frequency. A larger tip radius and a smaller 
tip-sample separation gives a more red-shifted resonance. 

Model calculations based on the physical picture de- 
scribed above gave results in rather good agreement with 
experiment, both in terms of absolute intensities and 
spectral properties. In these calculations, the compli- 
cated, extended geometrical shape of the STM tip was 
approximated by a sphere, and the electromagnetic field 
near .the .model tip was calculated in the non-retarded 
limit.E!rE3 The justification for the approximations is that 
the region between the tip and sample is most important 
for the effect. There the tip shape is well approximated 
by a sphere, and all relevant distances are small compared 
with the wavelength of the emitted photons. 

However, in view of future developments, it is rather 
useful to know how good these approximations are. In 
this paper, we focus on the effects of retardation. The cal- 
culations are built on exact diffraction theory, and leads 
to equations that can be solved numerically with good 
accuracy without too much effort. 

A number of calculations have addressed related 
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problems One example of .this., is scanning near- 
field optical microscopy (SNOM),HEE3 where one stud- 
ies photon tunneling from a transparent sample to a tip. 
The tip-sample distance may be several hundred A in this 
case, so it becomes necessary to take retardation into ac- 
count. But the calculations can often be carriedpOjiit us- 
ing various discretization schemes in real space,E3 since 
the relevant length scales differ by at most an order of 
magnitude. It deserves to be pointed out that the calcu- 
lational scheme developed in this paper could, with some 
extensions, also be applied to SNOM problems, and may 
then be faster than the prevailing methods. However, 
a prerequisite for this is that the tip-sample geometry 
possesses some (preferably cylindrical) symmetry. 

In anotheiJully retarded calculation, Takemori, Inoue, 
and OhtakaO (TIO) studied the field around a sphere in 
front of a plane surface with the aim of investigating elec- 
tromagnetiCj-effects in surface-enhanced Raman scatter- 
ing (SERS).Ej Finally, Madrazo, Nieto-Vesperinas, and 
Garciac3 (MNVG) have considered light emission from an 
STM in a calculation accounting for retardation. They 
used an essentially two-dimensional (2D) model geome- 
try, where the tip was represented by a horizontal cylin- 
der and the sample surface was artificially corrugated. 
We will compare our results with those of MNVG in as 
far as this is possible (they presented results for one single 
photon energy) and we will see that the different model 
geometries lead to quite different results. 

The inclusion of retardation effects does not change the 
basic picture of how light is emitted from an STM, how- 
ever, at a quantitative level, there are changes to some 
of the calculated results. As a general observation, retar- 
dation effects become increasingly important when the 
tip and/or sample materials are good conductors (such 
as silver). Then the energy dissipation in the near-field 
zone is rather small, and a localized mode will instead be 
damped because some energy propagates away from the 
near-zone. This is of course allowed for in a retarded, 
but not in a non-retarded, calculation. Thus for a Ag 
sample scanned by a Ir or W tip, retardation leads to 
an additional red-shift of the resonance frequency, and 
this improves the agreement between theory and exper- 
iment. For Au and Cu samples, on the other hand, the 
resonance frequency is essentially unchanged, and there 
are only relatively minor changes of the light emission 
intensity as a result of retardation effects. 

The rest of the paper is organized as follows: In Sec. ||, 
we develop the formalism needed to carry out a retarded 
calculation of the electromagnetic field that results from 
inelastic tunneling. We derive expressions for the radi- 
ated power, as well as the power dissipated in the sample 
and tip. In Sec. [II, we establish the connection between 
the completely classical current we use in the electrody- 
namic calculations, and the transition matrix elements 
associated with inelastic tunneling of an electron from the 
tip to the sample or vice versa. The numerical results of 
our calculation are presented in Sec. IV, and compared 
with experimental results. Finally, our conclusions are 



presented in Sec. |y|. 



II. THE ELECTROMAGNETIC FIELD 
CALCULATION 

A. General considerations 

We will calculate the electromagnetic field result- 
ing from excitation by a classical current distribution 
j e (r)e~ lult , localized between the tip and the sample. As 
in our previous, non-retarded calculations, we model the 
tip by a sphere with radius R that is centered at the ori- 
gin of our coordinate system, see Fig. [j]. The sample fills 
the half-space z < zq, where zq = —(d + R), d being the 
tip-sample separation. 

The optical properties of the tip and sample materials 
are described by local dielectric functions eii^) (sample) 
and €2^) (tip) determined from experiment H3 While this 
is an approximation, it is a reasonable one; we will discuss 
its limitations in some more detail in Sec. IV D. 



We write the source current distribution as 



j e (r) = zF(p)S(z - z )C. 



(1) 



The source is in other words a cylindrically symmetric 
dipole layer right at the sample surface. Most of the time 
we will take the radial form factor F(p) to be a Gaussian 



F( P ) 



= e~ p /a 



(2) 



with a ~5 A (as long as the current is concentrated to a 
region near the symmetry axis the exact shape of F(p) is 
unimportant). Even if we treat j e as a classical current, 
its strength C, must eventually be determined from a 
quantum- mechanical calculation, which we defer to Sec. 
III. 

Thanks to the cylindrical symmetry of both the cur- 
rent distribution and the tip and the sample, the mag- 
netic flux density is always directed along the azimuthal 
unit vector 0, and its magnitude only depends on z and 
the distance p, to the cylindrical symmetry axis, 



B(v) = B(p,z)cf>. 



(3) 



The accompanying E field is also cylindrically symmet- 
ric. In the region between the tip and the sample, it is 
mainly directed along z and thus strongly coupled to j e . 
Applying Maxwell's equations, we find that B must sat- 
isfy the vector Hclmholtz equation (we use SI units and 
po is the vacuum permeability) 



V 2 B + k 2 B = -noV x j e 



(4) 



in each of the three different regions of space. In vacuum 
feg = k 2 = cu 2 jc? . in the sample k\ — ei(uj)k 2 , and in the 
sphere fcf = £2(u)k 2 . Inside the sphere the solution to 
Eq. (||) can be written 
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ji{k 2 r) , ( 
>l^ a ' 777. u\ P i( cos9 )' 



1=1 



Ji(k 2 R) 



(5) 



where r, 6*, and are the usual spherical coordinates, ji is 
a spherical Bessel function, is an associated Legendre 
function, and a/ is a set of coefficients that remain to be 
determined. In the sample we write B as a Fourier-Bessel 
transform, 



(G) 



B = <f> I ndn B(k)Ji(kp) 



x exp 



-i\ jk\ - k 2 {z - zo) 



where J\ is an ordinary Bessel function, while B(k) is 
an as yet unknown function. Equation (ra) yields a finite 
field at the center of the sphere, Eq. (|^) gives outgoing, 
damped waves in the sample, and both expressions sat- 
isfy Eq. (Q). Next we must determine ai and B(k). 



B. The solution 



To calculate the field in the vacuum region, we use the 
vector equivalent of the Kirchoff integraH 



B(r) = B (e) (r) 



ik 



dS'[{n' x B(r')) x V'G G{n' x E(r'))]. (7) 

c 



Here, the first term B*- 6 -' solves Eq. (J4|) when neither the 
sample nor the tip is present, i.e., it yields the direct con- 
tribution to B from the source. G stands for the Green's 
function to the scalar Hclmholtz equation in vacuum, 



G(r,r') 



e ifc|r-r'| 
47r|r - r'l 



(8) 



The integration dS' in Eq. (|t]) runs over two surfaces: a 
plane just outside the sample and a spherical surface just 
outside the model tip. The normal vector n' points in the 
direction from the sphere or half-space into the vacuum 
region. The fields B(r') and E(r') are the exact fields on 
the two surfaces, so Eq. (|t]) is an integral equation. 

Before continuing, let us point out that even if the sur- 
face integral in Eq. ([?]) contains the full B and E fields 
on the surfaces, the combination appearing in the in- 
tegral guarantees that the final result only contains re- 
flected waves. If one for example wants to calculate the 
field from a certain source distribution in vacuum, but 
nevertheless introduces a surface somewhere outside the 
source, the integral does not at all contribute to B(r). In 
the present case, the integral over the spherical surface 
vanishes if e 2 = 1, and the integral over the flat surface 
vanishes if t\ = 1. 

Since Eq. (0) yields the field B at any point outside 
the tip and sample, it can be used to find expressions 



determining a; and B(k). Just outside the surfaces over 
which we integrate, we can (i) either calculate the field 
from Eq. (0) , or (ii) we can use Eq. (||) or Eq. (Q) because 
B is continuous across the interfaces. The so obtained, 
different expressions for B must be equal. Expanding the 
result outside the sphere in terms of associated Legendre 
functions yields 

/>oo 

ai = af^ + siai + / ndn (9) 
Jo 

where a[ e ' comes directly from the source, the term s/a/ 
comes from the integration over the sphere, and the 
last term originates from the sample surface integration. 
Without the sample, the above equation has the solu- 
tion a; = a, /(l — s/); thus, plasmon resonances occur 
on the isolated sphere when s; = 1. The corresponding 
equation, valid at the sample surface, is 



B(k) = B^(k) + S{k)B(k) + aigi(n). 



(10) 



i=i 



Removing the sphere, B(n) = B^(k)/(1 — S(k)), so the 
surface plasmon resonances occur when S(k) = 1. The 
explicit expressions for B^(k), //(k), gi(n), s;, S(n), 

of , and af are given below and in the Appendix. 

To solve for a; and B(k), we first note that Eq. ( fTo| ) 
has the formal solution 



B(k) 



1 - S(k) 



(11) 



Inserted into Eq. (Jj|) this yields a system of linear equa- 
tions determining the coefficients a;, 



where 



l'=l 



M w a v 



n ^ 1-S( K ) ' 



(12) 



(13) 



The terms on the right hand side are the driving forces 
and depend on the source: a, is the direct contribution 

from the source, while af^ results from the source fields 
reflected once off the sample surface. The two coefficients 
can be added together to give 



(e) , (p) 

a] + a\> 



00 h(n)B^(K) 
„ ^ S(k)[1-S(k)Y 



(14) 



B^ (k) is the Fourier-Bessel transform of the solution to 
Eq. (J|) at z — zq and with t\ — 1 and e 2 — 1- Using the 
expression in Eq. (|Al|) for the Green's function, we get 



2Vk 2 - k 2 Jo 
i[i^C na 2 



pdpF'ip)^) 



Vk 2 



i 2 4 



(15) 
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where the second line is valid if F(p) is given by Eq. (||). 

To summarize, a numerical calculation of the elec- 
tromagnetic field begins by evaluating si, Mu>, and 
from the formulas given here and in the Ap- 
pendix. Then the coefficients ai can be found by a matrix 
inversion in Eq. (p"2|). Having found a;, we get B(k) from 
Eq. (|l|), and subsequently B(r) can be calculated any- 
where in space. These steps must of course be repeated 
for each photon frequency u>. 

To carry out the matrix inversion in practice, a trun- 
cation in I at Z max must be introduced. With the geom- 
etry parameters that we consider most of the time (d ~ 
5A, and R < 500 A), l max ~ 50 is appropriate in order 
to get converged results. / max is determined by the size 
of the resonant cavity between the tip and the sample. 
The interface-plasmon mode is to a large extent confined 
within the region where the distance between tip and 
sample is less than twice the smallest tip-sample separa- 
tion d (cf. Ref. p8[ ). For a spherical tip this region has a 
radius y/2d~R (if d<R), i.e., with d = 5 A and R = 300 
A we get a "plasmon radius" of w 55 A. Since ttR w 950 
A, it is clear that even a crude description of the inter- 
face plasmon would require Z max ~ 20, but to capture its 
structure a larger number of basis functions is needed. 

The computation of the radiated power for one fre- 
quency takes less than 5 seconds on an ordinary work 
station. Most of this time is used for the k integration 
determining Mw in Eq. Jl3|). 



C. Emitted power 

Once the coefficients a; have been determined, we can 
evaluate the radiation power as well as the power dissi- 
pated into the sample and tip. In vacuum, far from the 
model tip, the B field has the asymptotic behavior 



ikr 



B 



'far" 



(16) 



-Bf ar is a sum of contributions from the source, the sample 
surface, and the spherical surface, 



5far= -[1 + S(k sin 6)}k cos 6B(k sin 9)e lk{d+R) cos 



^ ar k 

1=1 



4tt 



hi(kR) 



(17) 



Knowing i?f ar , we can calculate the radiated (differential) 
power per unit solid angle 



dPrad _ 1 3 , R , 2 



(18) 



By evaluating the Poynting vector on the surface of 
the model tip and then taking the time average, we find 
that the power dissipated into the sphere is given by 



P sph = ttR 



2 c 3 e 



E 

i=i 



21(1 + 1). a 



21 + 1 



aA Im 



1(1 + * 

C2 K R 2 ]i(k 2 R) 



(19) 



A similar calculation yields the power absorbed in the 
sample 



3 z 100 
p _ - C C ° / DC„M2 

J samp l / 

fc JO 



Kd K \B( K )\ 2 Re 



ei 



(20) 



D. Comparison with other methods 

Takemori, Inoue, and OhtakaEl studied scattering of 
a plane wave impinging on a sphere-plane system. This 
problem can be studied also in the present framework if 
B^ 6 -* is taken to describe the incoming wave. Of course, 
the solution would not in general be restricted to cylin- 
drically symmetric electric multipole modes as is the case 
here, but thanks to the cylindrical symmetry there would 
not be any mixing of modes with different m. The nu- 
merical solution would therefore be quite feasible. Note 
also that thanks to the reciprocity theorem,E] already the 
present calculation can yield results for the z component 
of the electric field induced by an incoming wave on the 
cylindrical symmetry axis below the tip. 

The conversion of plane waves to spherical waves and 
vice versa is central to both our solution and that of TIO. 
However, the conversion is dealt with in different ways. 
TIO essentially sum up the contributions from repeated 
scattering, back and forth, between the sphere and the 
plane. These repeated scattering events are implicitly in- 
cluded in the Kirchoff integral formulation. Thus, we feel 
that, at least once the field has been written in the form 
of Eq. (0), the present solution is conceptually simpler. 

Our calculation and-,thc one by Madrazo, Nieto- 
Vesperinas, and GarciaEElo have similar starting points; 
in both cases the electromagnetic field is expressed in 
terms of integrals over the surfaces where the relative 
dielectric function changes in a discontinuous way. How- 
ever, unlike us, MNVG choose to solve the resulting inte- 
gral equations in real space. While MNVG treat a situa- 
tion where the sample surface is periodically corrugated, 
it seems like their method would require quite intensive 
numerical calculations even for the case of a flat sample. 



III. THE SOURCE 

In order to perform the explicit calculations, we need 
to specify the sources. The tunnel current of an STM is 
typically concentrated to a small area; most of the cur- 
rent goes through the very last atom at the end of the 
tip. Also the inelastic transitions, from a state in one of 
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the electrodes to a lower-energy state in the other elec- 
trode, that drive the light emission are concentrated to 
the same small region of space. 

We wish to relate the classical ac current j e in Eq. (|l|) 
to quantum-mechanical transition matrix elements. Let 
us first look at one particular electronic transition from 
the. state \i) to |/). The corresponding classical current 
is§ 



j e W(r) = </|2j(r)|i>. 



(21) 



The resulting current distribution is of course not iden- 
tical to that of Eqs. (Q) and (0), however, the current is 
mainly directed along z and it is concentrated to a small 
region near the end of the tip. As a consequence, the two 
current distributions generate essentially the same fields 
provided their dipole moments 



d A rzp e {v) and / d A rzpi fi \r) 



are equal. Here the charge distributions p e and pi are 
connected to j e and j e through continuity equations 
such as 



(22) 



In order for the dipole moments to be equal, we must 
have 

- — ^ / pdpF(p) = — / d 3 r(f\2j z (r)\i). (23) 
iuj J Q iuj J 

Here C/i is the value C would take if (|i) — * |/)) was the 
only inelastic transition contributing to photon emission. 

In the experiment, light emission is caused by many 
different transitions at a whole range of frequencies. 
Therefore the factor C entering Eq. (Q) should be an 
incoherent sum of the different Cfi] moreover, we should 
calculate the emitted power per unit photon energy rather 
than just the emitted power. All the expressions for radi- 
ated or dissipated power found from (Eqs. (plf), ([l9|), and 
(f2p|)) are proportional to |C| 2 . We need to replace |C| 2 
by a new quantity, |C| 2 (w), in order to get expressions for 
the radiated or dissipated power per unit photon energy. 
For the transitions around photon energy Huj, we have 



\C\ 2 (u)=J2\CM 2 S(E l -E f -hco). 

fi 

With the aid of Eq. ( p3| ) , this can be rewritten as 

\(mM 2 



(24) 



|C| 2 H = E 



fi 



^\f™pdpF(p)\ 



: 8{E t -E f - Twj). (25) 



Now, for example, the radiated power per unit solid angle 
and photon energy is 



dP rad \C\ 2 (u;) 1 3 a 

-C Co -Bfar 



d{huj)dVL |C| 2 2 



(26) 



if the calculation of _Bf ar is still based on a source de- 
scribed by Eq. (|l|). 

In the actual calculation of |C| 2 (w), we have used free- 
electron models for both the tip and sample; more details 
can be found in Ref. ^2|. The resulting \C\ 2 is rather fea- 
tureless; the lj dependence is roughly 



KTM « i 



huj 

eVbim 



(27) 



where Vbias is the bias voltage. We have not tried to im- 
prove on these earlier calculations, for the effects-we are 
mainly interested in are of electrodynamic origin.Ea 

We should point out that we restricted the evaluation 
of the transition matrix elements to the vacuum (barrier) 
region. There are several good reasons for doing that: 
Since, in the frequency range where the photon emission 
is most intense, |etip| and |e S ampio| are considerably larger 
than 1, the electric field that couples to the tunnel cur- 
rent is stronger in the vacuum region. Moreover, the 
electron wave functions oscillate inside the sample and 
tip, whereas they do not change sign in the barrier re- 
gion. This means that the vacuum region should give the 
largest contributions to the matrix elements. 

Persson and Baratofli! studied this issue rather thor- 
oughly by comparing the contributions to photon emis- 
sion from inelastic tunneling and hot-electron decay. 
Within their model, in which electrons are tunneling 
into a spherical metal particle, inelastic tunneling pro- 
cesses (light emission while the electron is in the bar- 
rier) are about 3 orders of magnitude more effective than 
hot-electron decay (light emission when the electron has 
reached the final-state electrode) in the light-emission 
process. 



IV. RESULTS 



We have calculated the electromagnetic field for a num- 
ber of different materials combinations and geometric pa- 
rameters. Below we present results for the differential 
power (i.e. the radiated power per unit solid angle and 
photon energy) and the total radiated or dissipated power 
per unit photon energy. 

When we present results for the differential power, we 
take the observation angle to be 9 = 57.3° (i.e. 1 rad). 
The angular distribution of radiation shows only slight 
variations with frequency, materials, etc. It is given by a 
dipolar radiation lobe that has been "turned" away from 
6=90° to have a maximum at 0=55-60°. Very similar 
angular distributions show up both in experiments (Fig. 
30 of Ref. E3) and other calculations (see Refs. p4| and 



5 



A. Au sample — W tip 

We begin by considering light emission from a Au 
sample probed by a W tip. Figure || presents results 
for the differential power from both retarded and non- 
retarded calculations. As for the spectral shape, these 
results should be directly comparable with experiment 
(typically photons are collected over a certain solid angle 
in the experiments). □ Here, retardation does not cause 
any qualitative changes of the spectra. For a tip radius 
i?=100 A, the retarded and non-retarded results cannot 
be distinguished from each other. When R is increased to 
300 A, there is a difference in terms of intensity between 
the retarded and non-retarded results, but the spectral 
shape and resonance frequency is nearly the same. These 
spectra are in quite good agreement with the experimen- 
tal ones.el 

When the sample consists of Au (and the same holds 
true for Cu) , the resonance frequency is not very sensitive 
to the tip-sample geometry. The reason is that, because 
of the onset of interband transitions, the real part of the 
dielectric function £a u increases very rapidly from 
at 2 eV to «-2 at 2.5 eV. Thus, even though changes in 
the geometry leads to changes in the resonance condition 
expressed in terms of a value for the sample dielectric 
function, the resonance frequency does not change much; 
here there is only a slight red-shift of the resonance when 
R increases froni 100 to 300 A. 

Umeno et alS3 observed light emission from a granular 
Au film, and found variations in the resonance frequency 
that they interpreted as the result of a varying tip curva- 
ture. That interpretation is in qualitative agreement with 
our results, but a quantitative comparison is not possible 
for at least two reasons: (i) In its present form, our the- 
ory does not account for the granularity of the film, (ii) 
The bias voltage was as low as 2 V in the experiment, so 
the resonance frequency was most likely limited by Vbias • 
It would be interesting to see the experiment repeated at 
a larger bias. 



B. Ag sample — Ir tip 

Next let us look at the radiation spectra from a silver 
sample scanned by an Ir tip. In Fig. we display results 
of calculations including retardation with three different 
tip radii (200, 300, and 400 A, respectively), and for com- 
parison also results from a non-retarded calculation with 
i?=300 A. Here, the retarded and non-retarded results 
still have the same basic features, but the inclusion of 
retardation effects leads to substantial changes at the 
quantitative level. The primary resonance frequency is 
red-shifted, the second peak (between 3 and 3.5 eV) re- 
duces to a shoulder, the peak height is reduced by almost 
a factor of 2, and the peak becomes broader. 

One important reason for the lowering of the resonance 
frequency and increased damping (width) of the reso- 



nance appears to be that the fields penetrate further into 
the sample and tip when retardation effects are included 
in the calculation. This increases the resistive losses suf- 
fered by the tip-induced plasmon mode. We would like 
to stress that radiation damping as such plays no role in 
broadening the resonance here; as we will see below, the 
radiated power is just a small fraction of the dissipated 
power. 

The changes brought about by the retardation effects 
improve the agreement between calculated and exper- 
imental spectra. The second peak, seen in the non- 
retarded calculation but not in experiment, is gone, and 
the extra red-shift found in the retarded calculation 
brings the remaining resonance closer to its experimental 
position (typically 2.5 eV).B At the same time, a word 
of caution is in place here. The calculated spectra are 
clearly quite sensitive to the tip-sample geometry. In ex- 
periments, the peak positions varied somewhat between 
spectra taken with different tips on silver samples.u How- 
ever, these variations were not as large as those seen in 
Fig. |§|, where the resonance frequency is lowered by «0.5 
eV when R is increased from 200 A to 400 A. Unfor- 
tunately, our limited knowledge about the actual shape 
of the STM tip makes it difficult to reach any definite 
conclusion about this issue. 

Further insight can be gained by studying how the to- 
tal radiated power, as well as the power lost to the tip 
and sample, vary with frequency. Such results are dis- 
played in Fig. H for the case of a Ag sample and an Ir 
tip with i?=300 A. The losses to the tip and sample are 
much larger than the radiated power (note that the P ra d 
data are multiplied by a factor of 100). From the figure, 
the total radiated power can be estimated to 

25 x 10 6 y x leV w 4pW, 

or 10 7 photons/sec. The tunnel current 10 nA corre- 
sponds to w 0.62 x 10 11 tunneling electrons per sec- 
ond. Thus, the quantum efficiency is about 2xl0 -4 pho- 
tons/tunneling electron. This is in reasonable agreement 
both with the results otour earlier calculationsclrtj and 
experimental estimates.Bo 

Here one should note that the electromagnetic response 
is not always the limiting factor for the quantum effi- 
ciency. To reach the values found above, the bias voltage 
must be 1 V, or so, larger than the resonance photon en- 
ergy. But if the bias voltage is just 2 V, and in addition 
a Au sample and a W tip is used, the quantum efficiency 
drops to 10~ 6 or less (this is typically what was found 
in the experiments in Refs. || and |l0| ). This estimate is 
found by noting that: (i) if the bias voltage is 2 V, in 
view of Eq. (E7|), less than 10 % of the area under the 
curves in Fig. ^remains, and (ii) comparing Figs. ^]and|3| 
one realizes that the quantum efficiency for a Au-W con- 
figuration is about an order of magnitude smaller than 
for a Ag-Ir configuration. 

Returning to Fig. |I] we see that the dissipated power 
is about 100 times as large as the radiated power; thus, 
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since about 2 out of 10 4 tunneling electrons take part in 
photon emission processes, it is clear that at most a few 
percent of all the electrons undergo some inelastic pro- 
cess while in the gap between the tip and sample. (The 
calculation of Ref. [2?], using a two-sphere model geome- 
try, showed that the inelastic current could account for 
up to 10 % of the total current.) In any case, the to- 
tal tunnel current is dominated by the elastic part and 
the channels opened by inelastic processes only account 
for a rather small correction to the current. Most of the 
energy supplied by the bias voltage is eventually dissi- 
pated in hot-electron decay processes well inside the tip 
or sample. 

These observations are important because our scheme 
for determining the tip-sample distance d that gives a 
certain tunnel current J at a certain bias vol tag e Vbi as is 
now justified. In that calculation (see Ref. g2|) we take 
J = /clastic- A priori, nothing guarantees that the elastic 
contributions to the tunnel current are the dominating 
ones; given a sufficiently strong coupling to some other 
degrees of freedom, inelastic tunneling could dominate. 

In the present case, the resistive losses into the tip (Ir) 
dominate for most frequencies; Ir is a rather bad conduc- 
tor. However, just above 3.5 eV the losses into the sample 
are larger. This is mainly due to surface plasmon emis- 
sion (the surface plasma frequency of silver is ss 3.7 eV). 
Then the tunneling electron excites a surface plasmon 
that can propagate quite far away from the tip-sample 
cavity, but eventually loses its energy due to dissipation 
into the sample. 

Let us finally see what happens when the tip and sam- 
ple materials are interchanged. Figure [| compares the 
emission rates for a Ag sample-Ir tip with that of the 
combination Ir sample and Ag tip. In this, calculation 
we used a "white" spectrum for the source.E3 As is seen, 
the magnitude of the emitted power can change by a fac- 
tor of 2-3 as a result of the interchange in our model. 
Considering the details of the spectra there are further 
changes. The spectrum obtained with a Ag model tip 
has two peaks, and the one at higher frequency has the 
largest magnitude. Clearly, this peak originates, at least 
in part, from the fact that the model tip is a sphere, 
since it appears not very far below the frequency where 
an isolated silver sphere has a resonance. 



C. Ag sample — Ag tip 

The effects of retardation become even more apparent 
when we consider a situation in which both the tip and 
the sample are made of silver. Experiments have been 
done with tungsten tips covered by silver .□ The obtained 
spectra typically had two resonances. However, since the 
exact "composition" of the tip in these experiments is 
unknown, we will not calculate absolute intensities or 
make any comparison with experiment, but rather look 
at trends and make relative comparisons. 



In Fig. we display results for the radiation ampli- 
fication factor, i.e. the actual differential power divided 
by the value it would take if the tip as well as the sam- 
ple was absent, comparing the non-retarded and retarded 
calculations. Silver is a very good conductor over the en- 
tire frequency range up to huj—3.5 eV, and tip-induced 
plasmons localized to the region around the tip-sample 
gap are very weakly damped, unless energy can be dissi- 
pated due to wave propagation. Therefore the spectrum 
resulting from a non-retarded calculation has a series of 
very sharp peaks. Once the finite speed of light is taken 
into account, these resonance peaks are broadened and 
their heights are reduced by about an order of magni- 
tude. Still the amplification (enhancement) is very high. 
For example, if the corresponding curve was plotted for a 
Ag sample and an Ir tip, the peak value would be «10 5 . 
The changes brought about by including retardation can 
be understood by essentially the same reasoning as we 
used before. 

In Fig. ^, we show how the differential power (caL. 
culated with a white power spectrum for the currcntcfj 
develops with increasing tip radius. Each spectrum has 
several, rather sharp peaks. The different resonances cor- 
respond to modes with different field distribution in the 
gap between the tip and the sample. The resonance with 
the lowest frequency is nodeless, the second peak comes 
from a mode that has one node, etc. As before, the reso- 
nances are red-shifted and broadened when R increases. 
As for the peak height, it first grows because a larger 
sphere functions as a larger "antenna" putting out more 
power. However, already at i?=300 A the peaks are lower 
and broader due to the increased damping. When this 
happens, the total radiated power remains approximately 
constant to begin with. But eventually also this trend is 
broken once the size of the sphere becomes comparable 
to A ( A/27r). Then different parts of the sphere cannot 
send out radiation coherently any longer. 

The differential power spectra are of course also sensi- 
tive to the tip-sample separation d. This is illustrated in 
Fig. |^. Keeping R fixed and increasing d, the resonances 
are shifted upwards in frequency due to a weaker cou- 
pling between the tip and sample. At the same time, the 
peak values decrease; with a larger tip-sample separation 
the field enhancement in the cavity becomes less effec- 
tive, and the total photon yield decreases monotonically 
with increasing d. If one just looks at the intensity at a 
certain fixed photon energy, however, it can happen that 
the light intensity first increases as the tip approaches the 
surface, but then decreases again. This is due to the shift 
of the resonances as the tip-sample distance is changed; 
for a certain d a resonance appears at the frequency in 
question. An example of this is seen around tku = 2.1 eV 
in Fig. ^. 

In Fig. ||, we have plotted the total radiated and dife 
sipated power calculated with a white power spectrumo 
for the current. There is a conspicuous peak in P sam p 
just above 3.5 eV, which is due to surface-plasmon exci- 
tation. However, over most of the rest of the frequency 



7 



range, the radiated power is comparable to P S am P and 
P sp h- Thus, in the present case, radiation damping re- 
ally gives a significant contribution to the peak width. 



D. Surface modifications to the dielectric response 

As mentioned earlier, using dielectric functions de- 
termined from optical experiments is an approximation 
(however, as we will see, a good one) in the present prob- 
lem. A more accurate dielectric function should include 
effects of extra damping due to surface scattering and 
electron- hole pair excitations. Here we will try to give 
an idea about the size of this damping using rather sim- 
ple estimates. 

To estimate the change in e due to surface scattering, 
let us assume that the dielectric function at least over 
a certain frequency range can be written on the Drude 
form 



e(cj) = 1 - (-f ) +i— f- 



(28) 



where uj p is the plasma frequency, set by the density of 
conduction electrons, and r is the electron mean free 
time. A simple estimate of r is just r = £/vf, where 
£ is a mean free path and vf the Fermi velocity. Extra 
scattering off the surface should primarily occur in the 
tip, because it has of course a different shape than the 
flat surfaces used in the optical measurements in which 
the local dielectric functions were determined. But the 
standard tip materials Ir and W already have dielectric 
functions with large imaginary parts so an extra damping 
mechanism there makes no big difference. Let us there- 
fore look at Ag, for which u> p « 9eU,e3 and vf ~ 1.4 x 10 6 
m/s. With I = 300 A (i.e. the tip radius we used most 
often) and Tiuj — 2.5 eV, we get an increase in Ime 



A(Im[e]) = (-?-) 



1 



2.5' 4.6 x lO 1 ^- 1 x 2.1 x 10~ 14 s 



0.13. 



(29) 



In Fig. we show how an extra contribution (here 
taken to be frequency-independent) to Ime affects the 
radiation spectrum for a Ag-sample-Ir-tip combination. 
Extra damping leads to a decrease in the differential 
power, but the changes are still relatively small for re- 
alistic modifications (i.e., eAg+iO.l and eA g +*0.3) of the 
dielectric function. Maybe the best justification for using 
optically measured dielectric functions is provided by a 
comparison between Figs. || and [l(]; changes in the tip 
geometry have larger effects than surface scattering. 

Next we estimate the damping rate that result from 
the fact that a plasmon, which is a coherent superposi- 
tion of electron-hole pair excitations, under certain cir- 
cumstances can decay into real electron-hole pairs. To 
discuss this in a simple way, we will use the so called d- 
parameter theory of electromagnetic surface response of 



a jellium Ji9~0 The frequency-dependent (and complex) 
function d_\_{u>) tells where the centroid of the induced 
screening charge is situated relative to the jellium edge. 
The order of magnitude of d± is set by kp 1 , i.e. d± ~ 1 A, 
The surface response function of the jellium, g{q\\ , w),E2l 
is modified (to lowest nrder in q\\d±) from its classical 
value (e-l)/(e+l) to@ 



e(w + 1 



e(w) + 1 



2qud±(Lo) 



(30) 



Then, provided that q\\ is large enough that retarda- 
tion effects can be neglected (qii w /c) and using 
e(uj) = 1 — ujp/oj 2 , the surface plasmon at an isolated 
jellium surface has the frequency 



w =w s [l — q\\d±_{(jj s )/2\ 



(31) 



where the classical surface plasma frequency io s — 
LOp/^/2. Thus, for large enough qu, the resonance fre- 
quency will have a shift away from u> s depending on 
Re d± , and the damping of the surface plasmon depends 
on Im d±. The analysis leading to Eq. ( |3l| ) can be gen- 
eralized to the situation where two jellium surfaces, a 
distance L apart, are facing each other. The resonance 
frequency for the low-frequency interface plasmon (from 
which the tip-induced mode illustrated in Fig. [l] can be 
derived) is given by 



-i\\ L 



l + e~ q W L <?||d_L(wip) 



1 - e~ q w L 



, (32) 



if the two jellia are identical. 

Let us apply this formula to the case of the tip- induced 
plasmon. Using the parameter values Hu>ip = 2.5 eV, 
qii = 7r/100 A -1 corresponding to a typical plasmon di- 
ameter (cf. Sec. II B), and L = 5 A, we find that the 
electron-hole pair excitations will give an additional full- 
width-half-maximum broadening 

FWHM c _ h pairs = -2Im[7kt>ip] w 0.15 eV. 

To get this value, we used 

Im d±(uj/up = 0.3) « 0.16A, 

as calculated by LiebschEl for r s — 3 appropriate for sil- 
ver. It is clear that the estimated FWHM is considerably 
smaller than the resonance peak widths found in Figs. || 
and [H]. 

V. CONCLUSIONS AND OUTLOOK 

In this paper, we have presented a fully retarded cal- 
culation of the light emission rate from a scanning tun- 
neling microscope probing a noble metal surface. The 
major aim was to investigate in what way, and to what 
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degree, the inclusion of retardation changed the results, 
compared with our earlier non-retarded calculations Ej'E3 

In general, the results of the previous calculations still 
remain valid from a qualitative and semi-quantitative 
point of view. Light emission is resonantly enhanced due 
to the formation of a tip-induced plasmon mode in the 
cavity formed between the tip and sample. The enhance- 
ment amounts to 4 to 5 orders of magnitude compared 
with light emission in vacuum, and the quantum effi- 
ciency is typically 10 -4 photons per tunneling electron. 
The resonant photon energy is «2.1 eV for Au samples 
and «2.0 eV for Cu samples, i.e., essentially the same 
results as were found before. In the case of Ag samples, 
taking retardation into account lowers the resonant pho- 
ton energy by a few tenths of eV. 

Our results can also be comparecLwith some earlier caL 
culations addressing either SERS,L3 or light emission.Ej 
Note that this comparison is possible thanks to the reci- 
procity theoremHj The present results give spectra that 
are very similar to the ones obtained in Ref. [33]. Also 
in terms of intensities, the two calculations give similar 
results. A direct comparison of the numbers is not pos- 
sible; here the relevant quantity is the square of the field 
enhancement at one point in space, while in SERS one 
calculates the fourth power of the field enhancement av- 
eraged over a certain surface (a model sphere in Ref. 33). 

The results of a previous retarded calculation by 
Madrazo, Nieto-Vesperinas, and Garciac3 addressing 
light emission from an STM differ from ours on a couple 
of points: (i) MNVG find lower enhancement factors, 10 3 
at most, while we get 10 4 (Au sample-W tip), 10 5 (Ag- 
Ir), or even 10 6 (Ag-Ag). With a Ag sample and a Ag 
tip, and exactly the same parameter values for photon fre- 
quency, observation angle, and tip radius, and with d =8 
A, we get an enhancement that is about 50 times larger 
than that found by MNVG. Moreover, MNVG state that 
a major part of the enhancement (approximately a fac- 
tor 100) originates from the artificial corrugation of the 
sample. An important reason for the lower enhancement 
appears to be their using a two-dimensional geometry; 
their model tip is a horizontal cylinder and not a sphere, 
so the concentration of the electric field to the "end" 
of the tip is much less accentuated in their calculation, 
(ii) MNVG also found that interchanging tip and sam- 
ple materials changed the photon yield by as much as a 
factor of 20. They only found a considerable enhance- 
ment when the sample material is a noble metal. In the 
present calculation, interchanging tip and sample materi- 
als changed the photon yield by at most a factor 2-3 (cf. 
Fig. H). Again, this shows that the electromagnetic tip- 
sample interaction only plays a minor role for the field 
enhancement within the model of MNVG, whereas it is 
paramount in our model. 

The present theory can be extended in a number of 
ways. If the sample is covered by an overlayer, or the 
tip is covered by a spherical shell of a different mate- 
rial, the calculations can be carried through in essen- 
tially the same way after a few modifications have been 



made. With some more effort it should also be possible to 
treat cases where the tip has a more complicated shape, 
like for instance, a spheroid. In principle also regularly 
corrugated or randomly rough sample surfaces could be 
treated within an extension of the present framework, but 
the amount of numerical computations would be consid- 
erably larger than now. Finally, a calculation that could 
really treat an extended tip (i.e. one that continues out 
towards a tip-holder) would be useful. This is, however, 
a more difficult problem than those mentioned before. 

In addition to this, the experiments by Alvarado and 
co-workers on the polarization properties of the emitted 
light and its dependence on the surface and tip proper- 
ties have raised some interesting questions. For example, 
when the direction of magnetization of a Co film is re- 
versed, the degree of polarization of the emitted light can 
change by up to 20 %£3 For a magneto-optic effect, this 
is a very large number, and the effect could be impor- 
tant if it turns out that it can be used to study surface 
magnetism with good spatial resolution. CHowever, in 
a similar experiment on iron, Pierce et al!Ei found es- 
sentially no light polarization due to magnetism.) The 
experimental results of Ref. [l7| indicate that tip-sample 
interactions play a role also for this effect. The emitted 
light was more polarized in the tunneling regime (low 
bias) than in the field-emission regime. Moreover, it is 
possible that the shape of the tip is important in these 
experiments. A very high degree of polarization can be 
reached using an-asymmetrically shaped tip on a non- 
magnetic surface.Ej 

In summary, we have shown that retardation does play 
a role in light emission from STM's, but that the changes 
are not too big. Thus, theoretical studies of new, more 
complicated effects like the light emission from magnetic 
surfaces can very well be carried out neglecting retarda- 
tion. Comparing our results with those of Ref. |2^, we 
conclude that it is essential to retain a three-dimensional 
tip geometry in order to get an enhancement of the light 
emission rate that is comparable to those found in exper- 
iment. 
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APPENDIX: 

In this appendix we give explicit expressions for the 
functions /;(«;), <7z(k), S(k), and s;; we also indicate how 
these expressions are derived. 

To this end it is useful to express the Green's function 
of Eq. (^J) in representations suitable for cylindrical and 
spherical coordinates. In cylindrical coordinates 



00 2 — S 

G(r, r')= i ^ — A ™' cos (m((f> — <f>')) 

m=0 



-in 



ndn 



10 — k 

and for spherical coordinates 



Jm(np)Jm(np') ^k 2 -K 2 \z-z'\ 



I 



G(r,r') = ik^2ji(kr < )h l (kr>)J2 

(I - m)\ 
X (l + m)\ 



(2-iy )(2Z + l) 

47T 



/=0 m=0 

P/"(cos0)P/™(cos6>')cos(m(</>- 0')), (A2) 



where r (r ) is the lesser (greater) of r and r' , and hi is 
a spherical Hankel function. In addition, the integral 
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(A3) 



will be used repeatedly. 

Let us begin by looking at the contributions to B in 
Eq. (Q) coming from the sample. From Eq. (^]), 



B(r') 



K,dKB(n)Jl(Kp') 



(A4) 



at the sample surface. To calculate the electric field just 
outside the sample surface, we first evaluate the E field 
inside the sample from 



E 



ic 



V x B, 



(A5) 



and then use the fact that the tangential component of 
E is continuous across the sample surface. This yields 



-in 1 x E(r')) = -« 
c ei 



ndnB(K)Ji(Kp')\/ k\ — k 



(A6) 



at the sample surface. Inserting Eqs. (A4), (A6), and 
( Al) into p) and carrying out the surface integral, we 
find that the sample surface contribution to B(r) every- 
where in the vacuum region can be written 
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B (s) (r) =4> KdKB(K)Ji(Kp)< 
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(A7) 



Comparing this with Eq. (p_0[), it is immediately clear 
that 



S(k) 



ei\/k 2 — k 2 



(A8) 



Furthermore, by calculating the overlap at the spheri- 
cal surface between B^(r), as given by Eq . (A7), and 
P/^cosfl), one finds, with the aid of Eq. (A3), the plane- 
to-sphere coupling 



21 + 1 -'-5(«) 



1(1 + 1) 

xe iVW^(d+R) p i (^/i-tf/k*} ji(kR). (A9) 



Note that if t\ — 1 both S(k) and fi(n) vanish. 

Next, we apply the same methods to calculate the con- 
tributions to B(r) from the spherical surface. This yields 



si= ikR hi(kR) ji(kR) 



x kR 



3l(kR) 



1 

1 

£2 



k 2 R 



W2R) 



ji(kR) 

for the sphere-sphere coupling, while 
i l R 2 



ji(k 2 R) 



(A10) 
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'ji(kR) " r R e 2 



K 2 V / 

MkR) (All) 



1 1 



1 , h ms) 

R ji(k 2 R) 



describes how the waves originating from the sphere 
propagates to the sample surface. In this case, £2 = 1 
causes si and gi(n) to vanish. 
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FIG. 1. Illustration of the tip-sample geometry (not to scale). In (a) the current j e exciting the system and the electric 
field associated with the tip-induced plasmon mode are shown. In (b) the model geometry and coordinate system employed in 
the calculations are shown. 

dP/d(huj)dn(10 6 s- 1 ) 




1.6 1.8 2 2.2 2.4 2.6 2.8 3 
Photon energy (eV) 

FIG. 2. Radiated power per unit solid angle and photon energy for the case of a gold sample and tungsten tip. The bias 
voltage is 3 V (negative sample) and the tunnel current is 10 nA. The corresponding tip-sample separation d was 5.71 A (#=100 
A) and 6.30 A (7?=300 A). The continuous curves show results of retarded calculations, while the dotted curves show results 
calculated without accounting for retardation. 
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FIG. 3. Radiated power per unit solid angle and photon energy with a silver sample and iridium tip. The bias voltage is 4 
V, the current is 10 nA, and the corresponding tip-sample separations are: d=6.09 A (for #=200 A), d=6.28 A (#=300 A), 
and d=6.41 A (#=400 A). As in Fig. 0, the continuous (dotted) curve shows results from a retarded (non-retarded) calculation. 
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FIG. 4. The total power per unit photon energy emitted as radiation, and dissipated into the tip and sample. The calculation 
takes retardation into account. The tip and sample materials are Ir and Ag, respectively, the tip radius #=300 A, and the 
rest of the parameter values are the same as in Fig. |^. To facilitate direct comparison, the radiation power data have been 
multiplied by a factor of 100. 
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FIG. 5. Effects of interchanging tip and sample materials. The two curves show the calculated differential power, given 
that the driving current has a white spectrumO 
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FIG. 6. Radiation amplification factor for an Ag sample and Ag tip with _R=300 A and d~5 A. The amplification factor is 
the ratio between the radiated differential power with tip and sample present, and the radiated power into vacuum given the 
same source distribution. The dotted curve gives the results from a non-retarded calculation. The peaks of this curve are about 
10 times higher and sharper than those of the continuous curve resulting from a retarded calculation. Note the logarithmic 
scale. 
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FIG. 7. Results for the differential power for the case of a Ag sample and Ag tip. A white power spectru was used for 
the current driving the light emission. The calculations were done taking retardation into account. The tip radii were 100, 200, 
300, and 400 A, respectively, as indicated next to the curves, and d=5 A. 
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FIG. 8. Differential power (calculated with a white spectrumB for the driving current) for a Ag sample and Ag tip with 
i?=300 A for varying tip-sample distance. d=5 A for the curve with highest peaks, d=7.5 A for the dotted curve, and d=10 A 
for the curve showing the lowest intensity. 
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FIG. 9. Total radiated, and dissipated power per unit photon energy for a Ag sample and Ag tip (R=300 A and d=5 A). 
The thin curve gives the total radiated power, the other curves showrfche dissipated power into the sample (thick curve) and 
tip (dotted curve). The results were calculated with a white spectrumCj for the current. 
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FIG. 10. Differential power (calculated with a white spectrumCJ for the driving current) for a Ag sample and Ir tip with 
i?=300 A and d=5 A. The different curves show how the spectrum develops as the imaginary part of the Ag dielectric function 
is increased: The thick curve gives the result obtained with £A g determined from optical experiments. The remaining curves 
were calculated using £A g + iO.l, £A g + i0.3, and eA g + il.O, respectively. 
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